Search for fingerprints of disoriented chiral condensates in cosmic ray showers 
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Although the generation of disoriented chiral condensates (DCCs), where the order parameter 
for chiral symmetry breaking is misaligned with respect to the vacuum direction in isospin state, is 
quite natural in the theory of strong interactions, they have so far eluded experiments in accelerators 
and cosmic rays. If DCCs are formed in high-energy nuclear collisions, the relevant outcome are 
very large event-by-event fluctuations in the neutral-to-charged pion fraction. In this note we search 
for fingerprints of DCC formation in observables of ultra-high energy cosmic ray showers. We 
present simulation results for the depth of the maximum (Xmax) and number of muons on the 
ground, evaluating their sensitivity to the neutral-to-charged pion fraction asymmetry produced in 
the primary interaction. 

PACS numbers: 



I. INTRODUCTION 

Almost forty years ago exotic, apparently hybrid and 
unexpected events, dubbed Centauros, were observed in 
cosmic ray (CR) experiments in emulsion chambers in 
Chacaltaya by Lattes and collaborators [I] . Those events 
were very different from what is commonly observed in 
CRs, exhibiting a large number of hadrons and a small 
number of electrons and gammas, which suggests the 
presence of very few rapid-gamma-decaying hadrons. So, 
a possible imbalance in the number of neutral to charged 
pions could be envisaged. The nature and reality of Cen- 
tauro events started a long debate, that includes the re- 
examination of the original emulsion chamber plates, and 
is still not resolved [2H5]. 

Nevertheless, Centauro events were certainly an exper- 
imental motivation for the development of the theory of 
disoriented chiral condensates (DCCs) that started in the 
early 1990s For a detailed review, see Ref. [TU] . 

Assuming that a given nuclear system could be heated 
above the critical (crossover) transition region for chi- 
ral symmetry restoration, i.e. for temperatures of the 
order of 180 — 190 MeV [TT], then quenched to low tem- 
peratures, the chiral condensate initially melted to zero 
could grow in any direction in isospin space. Besides 
the vacuum (stable) direction, it could build up as in a 
metastable, misaligned pseudo-vacuum state, and later 
decay to the true, chirally broken vacuum. The fact that 
DCCs could be formed in high-energy heavy ion collisions 
stimulated several theoretical advances and experimental 
searches [10]. Most likely the temperatures achieved in 
current heavy ion experiments are high enough to pro- 
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duce an approximately chirally symmetric quark-gluon 
plasma, and the following rapid expansion can cool the 
system back to the vacuum [12], so that the dynamics of 
chiral symmetry restoration and breakdown can be de- 
scribed in a quench scenario [S], so that the evolution of 
the order parameter can be much affected by an explo- 
sive behavior that naturally leads to large fluctuations 
and inhomogeneities [T^HTT] . 

Since, by assumption, the order parameter for chiral 
symmetry breaking, i.e. the chiral condensate, is mis- 
aligned with respect to the vacuum direction (associated 
with the cr-direction in effective models for strong inter- 
actions) in a DCC, this would be a natural candidate to 
explain the excessive production of hadrons unaccompa- 
nied by electrons and photons, suggesting the suppression 
of neutral pions with respect to charged pions. 

Regardless of the outcome of the debate on the nature 
of Centauro events, DCC formation seems to be a quite 
natural phenomenon in the theory of strong interactions. 
However, given its symmetric nature (in isospin space), it 
should be washed out by standard event averaging meth- 
ods. So far, there has been no evidence from colliders 
or CR experiments. Motivated by the possibility of at- 
taining much higher statistics in current ultra-high en- 
ergy cosmic ray (UHECR) experiments than in the past, 
so that an event-by-event analysis for very high-energy 
collisions can in principle be performed, we consider pos- 
sible signatures of DCC production in CR air showers. 
If DCCs are formed in high-energy nuclear collisions in 
the atmosphere, the relevant outcome from the primary 
collision are very large event-by-event fluctuations in the 
neutral-to-charged pion fraction, and this could affect the 
nature of the subsequent atmospheric shower. Very pre- 
liminary, yet encouraging results were presented in Ref. 

m- 

In this paper we search for fingerprints of DCC forma- 
tion in two different observables of UHECR showers. We 
present simulation results for the depth of the maximum 
(Xmax) and number of muons on the ground, evaluat- 
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ing their sensitivity to the neutral-to-charged pion frac- 
tion asymmetry produced in the primary interaction. To 
model the effect from the presence of a DCC, we sim- 
ply modify the neutral-to-charged pion fraction, assum- 
ing that the system follows the same kinematics, as will 
be detailed below. Although this is certainly a very crude 
description of the dynamics of the primary collision, we 
believe it captures the essential features that have to be 
tested in order to verify the feasibility of detecting DCCs 
in UHECR showers. 

This paper is organized as follows. In Section II we 
briefly review some characteristics of DCCs, especially 
the Baked-Alaska scenario and the inverse square root 
distribution of the neutral pion fraction. In Section III 
the method for the simulation is presented. We use COR- 
SIKA (19j , a program for detailed simulation of extensive 
air showers initiated by high-energy cosmic ray particles. 
In Section IV we show and discuss our results. Section 
V contains our conclusions. 



II. DCC FEATURES AND THE NEUTRAL 
PION FRACTION DISTRIBUTION 

It is widely believed that for high enough energy densi- 
ties, achieved e.g. by increasing dramatically the temper- 
ature, strong interactions becomes approximately chiral 
(it would be an exact symmetry only if current quarks 
were strictly massless), so that the chiral condensate, 
which is the order parameter for that transition, essen- 
tially vanishes. On the other hand, for low temperatures 
the chiral condensate acquires a non- vanishing value and 
breaks spontaneously the chiral symmetry of (massless) 
QCD ,20.. 

In a given model, one can construct an effective po- 
tential for the chiral condensate degrees of freedom and 
study the mechanism of chiral symmetry restoration and 
breakdown. If we restrict our analysis to two flavors of 
light quarks, up and down, that can be easily accom- 
plished by the linear sigma model coupled to quarks [H] . 
In that case, the effective degrees of freedom are pions, 
tt", and the sigma, a. In the high-temperature limit all 
field expectation values vanish, whereas in the vacuum 
one has (tt") — and (cr) = /^r, where is the pion 
decay constant. 

The physical picture we have in mind is a very high- 
energy heavy ion collision that will create a hot quark- 
gluon plasma where chiral symmetry is approximately re- 
stored. As the plasma is quenched to low temperatures 
by expansion, the system will evolve to the vacuum emit- 
ting a large number of pions. However, the evolution can 
proceed along many different paths in chiral space before 
it finally reaches the true vacuum, i.e. it can "roll" into 
different directions, and the ratio of neutral to charged pi- 
ons produced depends strongly on the chosen metastable 
state in each event. In other words, the misalignment of 
the vacuum is reflected in the distribution of produced 
pions, generating a coherent state and an anomaly in the 



ratio of charged to neutral pions. This effect will be, of 
course, washed out by event averages. 

A more intuitive physical picture for the formation 
of DCCs, the Baked-Alaska scenario, was proposed by 
Bjorken and collaborators ^[^Hl^^. Consider a high- 
multiplicity collision at high energy. Before hadroniza- 
tion, most of the energy released at the collision point is 
carried away by the primary partons at nearly the speed 
of light. This hot and thin expanding shell isolates the 
relatively cold interior from the outer vacuum. As the 
state evolves, the interior cools down and what was be- 
fore a global minimum becomes a local maximum. The 
symmetry is spontaneously broken and one of the pseudo- 
vacua should be chosen, with a slight preference for the 
true vacuum. However, if the lifetime of the shell is short 
enough, the quark condensate in the interior might be ro- 
tated from its ordinary direction, since it costs relatively 
little energy, e - f^ml « 20 MeV/fm^ [TO]. When the 
hot shell hadronizes, the reorientation induced by the ex- 
ternal contact is reflected in the produced pions, which 
will be coherent and present fluctuations in an event-by- 
event analysis. 

Assuming that the correlated region is large enough 
to be described semiclassically, one can use the linear 
sigma model with explicit symmetry breaking for the de- 
scription of the dynamics |22j. Quantitatively, one can 
represent the shell of this fireball as the source of these 
excitations of pions: 

{D+ml)7T'^{x) = r{x), (1) 

where a is the isospin index. A spherically expanding 
shell is represented by 

r{x) = r{t)Q{t-to)S{t-r), (2) 

with initial radius — to, where to is the time where the 
expansion starts. After the hadronization, the currents 
J"' (t) vanish and the fields 7r° decay towards the vacuum 
into freely propagating pions. So, the pion emission is 
characterized by the state 

|J> = AAexp ' (3) 

where A/" is a normalization factor, and the sum is over 
isospin directions a = 1, 2, 3. The creation operator of a 
pion with momentum k and isospin component a is cl^{k) 
and Ja{k) is the 4-dimensional Fourier transform of the 

source J°(x) at fc*^ = \Jk^ + rn^ . 

The number of pions follows a Poisson distribution 
with average 

n'\k) = {J\cl{k)cAk)\J)^\J''{k)\\ (4) 
and the number of pions produced per unit phase space 
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is approximately given by |10j 



The neutral pion fraction distribution, Eq. (10), is 



(J) 



(Pk 



(27r)3 



(5) 



Statistically, one expects that the magnitude and the 
chiral orientation of the source J° {k) will fluctuate event 
by event for each mode. Let us assume that the vacuum 
orientation is tilted into one of the pion directions, i.e.: 



{a) = /tt cos ( 



(6) 



so that all relevant modes of chiral condensate point in 
the same direction u in isospin space: 



Joccik) = J{k)u. 



(7) 



In generic models of production the neutral pion fraction, 
defined as 



(8) 



is a binomial distribution with average / = 1/3. In this 
way, in the limit of large numbers N^^a , the probability of 
all pions being charged is very small. However, if pions 
are the product of a DCC decay, this probability is not 
negligible. In fact, if there is no privileged isospin direc- 
tion, the vector ii can point in any direction within the 
unit sphere. Then: 



focc = cos^ e , 



(9) 



where 9 is the angle between the unit vector u and the 
TT^ direction. So, one finds the following well-known dis- 
tribution for the neutral pion fraction [7j |8] : 



dP w d(cos 6) = ;-^d(cos2 9) '^^ 



2 cos 6 



2Vf 



(10) 



The probability of less than 10% of pions be neutral, for 
instance, is ~ 30%. 



III. DCC SIMULATIONS 

The conditions of a high temperature initial state fol- 
lowed by a rapid cooling stage are both possible to hap- 
pen in heavy ion collisions at the very high energies like 
those accessible at RHIC and to be reached at the LHC, 
as well as in UHECR collisions in the top of the atmo- 
sphere. The large aperture of a detector like that of the 
Pierre Auger Observatory [24 (combining shower sam- 
pling at ground level and longitudinal shower profile re- 
construction) has been providing high quality data and 
unprecedented statistics in the field [25H27] . Therefore, 
even though the formation of a DCC is probably rare, 
we believe it is worth studying the implications of such 
events to the physics of showers generated by UHECR. 



at the basis of our strategy to search for DCC finger- 
prints in UHECR showers. Therefore, any investigation 
to measure the impact caused by the presence of a DCC 
should assess the sensitivity of a given observable with 
respect to the neutral pion fraction produced in the pri- 
mary interaction. This fraction determines the initial 
distribution of particles between the electromagnetic and 
hadronic components of the showers. In this paper, we 
consider two observables which are usually measured by 
UHECR detectors: the slant depth in the atmosphere 
(defined as the integral of the atmosphere density along 
the shower axis^ / pdx and expressed in units of g/cm^) 
at which the shower reaches its maximum development, 
^max J a-nd the number of muons on the ground iV^ . It is 
known that the parameter Xmax is affected by the first 
interaction cross-section and its associated multiplicity 
and inelasticity [5 . 

If DCCs really exist, the conditions for them to be 
produced should include not only high energy densities, 
but the regions where such densities are achieved should 
not be small as well, since DCCs are considered "macro- 
scopic" space-time regions where the chiral parameter 
is not oriented in the same direction as the vacuum. 
With those requirements in mind, we have chosen to 
work with Fe initiated showers at 10^^ eV. Then, central 
collisions are privileged over peripheral ones by select- 
ing events with a large number of participating nucleons 
[Npart > 40) which, in turn, should lead to a higher mul- 
tiplicity in the first interaction. 

For all simulations presented in this work, we have 
used CORSIKA 6.617 [19 with the interaction models 
Sibyll 2.1 28J and GHEISHA 2002d 29J, for high and 
low energy processes, respectively. The DCC-like shower 
simulation chain is as follows: Large Npart collisions are 
selected and their first secondaries (mostly pions and 
kaons) separated; after, some of the neutral pions in this 
sample are converted into charged ones; the resulting par- 
ticle list is then inserted back into CORSIKA to proceed 
with usual cascade development through the atmosphere. 
Such a procedure was performed for several tt'^ fractions 
and 2 different zenith angles. The first interaction slant 
depth (ATo) distribution for DCC-like showers will there- 
fore be the same as for normal Fe initiated showers in 
central collisions. This is a valid assumption, since the 
DCC formation process takes place during a subsequent 
cooling stage of the initial hot plasma, with the first in- 
teraction cross section being the same as in a Fe-nucleus 
collision. We believe that even though the simulation ap- 
proach adopted here is simplified, the essential features 
of the process are being taken into account. 



^ For brevity, from now on we should refer to the slant depth as 
simply depth. However, the reader should be aware that in the 
literature the term depth may be used to refer to the integral 
of the atmosphere density along the vortical and not along the 
shower axis. 
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For comparison, proton initiated showers as well as 
normal Fe initiated ones were also generated. For the 
normal Fe case, we have produced both a sample rich in 
central collisions and another sample with all the central- 
ities. For each shower we extract the value of Xmax and 
the number of muons on the ground. 



IV. DISCUSSION OF RESULTS 

From now on we shall identify our DCC-like Fe initi- 
ated showers by Fe+DCC. Fig. [T] shows the X,nax distri- 
bution of a vertical shower (6 = 0) corresponding to an 
extremely asymmetric situation (/ = 0), that is, where 
all the initially produced 7r°'s are converted into tt^/tt^. 



For the distribution of Eq. ( 10 ) , sharply peaked at / = 0, 
the probability for less than 1% of tt^'s being produced 
is 10%. Four types of showers: normal Fe (central col- 
lision: Fe-Central), Fe (all centralities) , Fe-|-DCC and 
proton initiated are shown. On can clearly see that both 
samples generated from a central collision have smaller 
than average Xmax-, since the higher the multiplicity in 
the first interaction the faster is the subsequent cascade 
development in the atmosphere. And there is essencially 
no difference between the Fc-Central and the Fe-DCC 
samples in terms of X^ax, which might indicate that 
the early stage where the 7r° population is depleted to- 
gether with the higher interaction probability due to the 
large multiplicity allow for a complete recovery of these 
particles in the subsequent interactions. Nonetheless, it 
is clear that one should look at low X^ax events when 
searching for DCCs signatures, and this property is in- 
dependent of the initial vr" fraction. 

Dislribution of Xmax I 
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centralites. This is to be expected, since the fluctua- 
tions in Xmax have basically two components: the the 
ones in the first interaction slant depth Xq (which, in 
turn, depends on the interaction cross-section) and those 
introduced by the cascade growth up to the maximum. 
The latter depends on the initial shower size (the mul- 
tiplicity), in average being smaller for high multiplicity 
events, i.e. those found in the central collisions. 
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Figure 2: Normalized distributions of the number of fj,'^ + 
fj,^ on the ground for vertical showers initiated by Fe (central 
collision), Fe-|~DCC, Fe and proton with E = 10^^ eV. 
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Figure 3: Depth of the maximum versus number of jj.'^ + 
fj,~ on the ground for vertical showers initiated by Fe (central 
collision), Fe-|-DCC, Fe and proton with E = 10^^ eV. 



Figure 1: Normalized Xmax distributions for vertical showers 
initiated by Fe (central collision), Fe+DCC, Fe and proton 
with E = 10^^ eV. 



Another feature which can be appreciated in Fig. [T] 
is the smaller X^ax fluctuations for Fe-DCC and Fe- 
Central as compared to those of proton and Fe with all 



The muon number distributions on the ground for the 
same showers described above are shown in Fig. [2j As 
expected, increasing the number of initial charged pions 
will lead to a corresponding larger density of muons on 
the ground as the products of the charged pions decay. 

As one combines both pieces of information {X^ax and 
Nfj_), one finds a good separation between the Fe-DCC 
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Figure 4: Depth of the maximum versus number of + fi on the ground for showers with 6 = 38 degrees initiated by Fe 
(central collision), Fe+DCC, Fe and proton with E = 10^^ eV for different 7r° fractions: /=0.0 (top left), 0.25 (top right), 0.50 
(middle left), 0.75 (middle right) and 1.0 (bottom left). Bottom right: merit factor (Eq. 11 1 as a function of /. 



sample and all the other populations, as can be seen 
from Fig. |3] When collisions with all the centralities are 
allowed, iron showers exhibit a positive correlation be- 
tween Xmax and Nfj_. The correlation is such that show- 
ers reaching maximum development at increasingly larger 
slant depths produce, in turn, a higher muon density at 
ground level, due to the decrease in the attenuation in 
the atmosphere from Xmax to the ground. 



muons is present^, as can be seen from Fig. [sj Such an 
inverted correlation can be understood remembering the 
additional correlation between Xmax and the multiplic- 
ity in the first interaction, with low Xmax events corre- 
sponding, in average, to high multiplicity. As we convert 
neutral pions to charged ones, high multiplicity (small 
Xmax) events show a larger muon density on the ground 
as compared to low multiplicity ones (large Xmax)- 



Whereas for Fe-Central showers one sees no correlation 
between Xmax and iV^, for the case of Fe-DCC showers, 
instead, a small anti-correlation of about 1.96 g/cir? / 10^ 



^ A fit to the binned X^ax X distribution provides (1.96 it 
0.13) X 10"'' g/cm^/muon. 
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Since the number of muons on the ground depend both 
on the TT^ fraction / and on the zenith angle, due to 
atmospheric attenuation, we decided to perform a more 
systematic study by varying those parameters. Applying 
a Linear Discriminant Analysis (LDA), similar to what 
has been done in [30] we have determined, as a function 
of /, a power of discrimination defined as 1 — /3, with the 
classification error /3, in the 2 populations case {A and 
B) given by: 



Nab + Nba 
Na + Nb 



(11) 



where Nab {Nba) are the number of events from popu- 
lation A (B) misclassified as B (A) and Na {Nb) is total 
number of events in population A (B). As training sam- 
ples we have used half of the simulated showers, with the 
other half being used for the determination of the power 
of discrimination. 

Fig. |4] (top, middle and bottom left) shows again the 
two-dimensional scatter plots in the Xmax x N^^ plane for 
the same 4 populations analyzed above and 5 values of 
7r° fraction / = 0,0.25,0.5,0.75,1 and 6* = 38 degrees. 
The merit factor as a function of / is shown at the bot- 
tom right plot of Fig. |4] One sees that for any value 
of /, proton showers can always be reasonably separated 
from Fe-DCC, with merit factors above 97% in the whole 
interval of /. For Fe showers, of course, as we approach 
the / = 1/3 fraction, the separation power gets weaker. 
Fe-Central showers are easier to separate from the Fe (all 
centralities) signals only at small values of /, where the 
excess of muons at ground is large enough to segregate 
the Fe-DCC population. Fe-DCCs showers in which the 
first interaction produces a large amount of tt" are very 
poor in muons and would also be easily recognized. How- 
ever, even if DCCs exist in nature, such large-/ events 
are very unlikely to happen taken into account the dis- 



tribution of Eq. ( 10 1 



A more robust estimation of the power of discrimina- 
tion would have to incorporate a priori information on 
the population frequencies for proton and Fe, as well as 
the DCC occurrence frequency. Nonetheless, since nei- 
ther the cosmic rays chemical composition at the highest 
energies {E > 10^^ eV) nor the full dynamics of DCC 
formation are known, flat priors on these variables are 
the best one can do so far. 



In this paper, using air shower simulations, we searched 
for fingerprints of DCCs formed in ultra-high-energy 
{E ^ 10^^ eV) central iron collisions in the upper at- 
mosphere. In particular, we studied the influence of the 
DCC formation on the observables X„iax and iV^, via the 
asymmetry in the neutral-to-charged pion ratio in the 
primary collision. For comparison, we considered also 
regular air showers generated by protons and iron. 

Since DCCs are expected to be formed in central colli- 
sions, which lead to smaller than average values of X^ax 
(a difference of about 40 g/cm^ as compared to iron, for 
instance), one should concentrate searches within this re- 
gion of shower depth. For the same reason, if DCC events 
are present they should yield small fluctuations in Xmax ■ 
However, based only on the behavior of this observable 
one can not distinguish between an iron shower and a 
central collision and one produced in the presence of a 
DCC. 

The formation of DCCs is associated with large event- 
by-event fluctuations in the ratio of neutral to charged 
pious, /. In particular, for / large or small as compared 
to 1/3, one should expect a sizable change in the muon 
density on the ground, especially for vertical showers. 
This fact was noticed in the preliminary study of Ref. 
|18) . In the extreme case of / = 0, where DCC events 
lead to muon-rich showers, we showed that this signature 
distinguishes between the cases of iron (even for central 
collisions) and a DCC event. For large /, one can also 
separate these two cases due to the large depletion of 
muons on the ground. Even in this case, the signature is 
not contaminated by proton events. 

For vertical showers, there is a clear anti-correlation 
between Xmax and N^. This comes about since there is 
a correlation between Xmax and the first interaction mul- 
tiplicity, as was discussed in the previous section. This 
behavior is not expected in a regular iron shower. In 
fact, for iron, due to atmospheric attenuation, one would 
expect a positive correlation. 

Even though the analysis presented here is very sim- 
plified, it has the advantage of providing a setup that is 
totally under control for simulations, and we believe it 
contains the essential ingredients of the phenomena con- 
sidered. Nevertheless, a more realistic study, that should 
contain a description of the dynamics of DCC formation, 
especially its dependence on energy, is certainly neces- 
sary. 



V. CONCLUSIONS 

The recent increase in quality and statistics of UHECR 
data brings the possibility to look for exotic phenomena 
as well as rare phenomena within the Standard Model. 
Among the latter, DCCs have been predicted about 
twenty years ago, yet never detected due to their elusive 
character. 
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